Strongly Non-Arrhenius Self-Interstitial Diffusion in Vanadium 
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We study diffusion of self-interstitial atoms (SIAs) in vanadium via molecular dynamics simu- 
lations. The (lll)-split interstitials are observed to diffuse one-dimensionally at low temperature, 
but rotate into other (111) directions as the temperature is increased. The SIA diffusion is highly 
non-Arrhenius. At T < 600 K, this behavior arises from temperature-dependent correlations. At 
T > 600 K, the Arrhenius expression for thermally activated diffusion breaks down when the mi- 
gration barriers become small compared to the thermal energy. This leads to Arrhenius diffusion 
kinetics at low T and diffusivity proportional to temperature at high T. 



PACS numbers: 66.30.Fq, 61.72.Bb, 61.82.Bg 

The creation and migration of self-interstitial atoms 
(SIAs) are critical for microstructural evolution of mate- 
rials in a variety of situations, such as in the high energy 
radiation environment of nuclear reactors 1] and in ion 
implantation 0]. Although SIA formation energies are 
much larger than typical thermal energies, they form in 
abundance during collision cascades induced by imping- 
ing energetic particles. SIAs in metals are typically very 
mobile (i.e., their migration barriers are relatively small) 
and hence play an important role in controlling the rates 
of many microstructural processes in such applications, 
in particular the phenomenon of void swelling. 

Since SIA properties and mobilities are very difficult to 
determine experimentally, one often employs computer 
simulations 0, 0, IE 0- For example, simulations of 
body centered cubic (bcc) iron (and several other bcc 
metals), have shown that SIAs preferentially lie along 
(110) orientations but rotate into (lll)-directions, where 
they can migrate easily using the crowdion configura- 
tion as transition state. Other simulations have sug- 
gested that SIA migration in vanadium is very similar 
to that in Fe 0, S S ^| ■ However, these empirical in- 
teratomic potential-based simulations are not consistent 
with recent first principles calculations that clearly show 
that the lowest energy SIA configuration in V is a (111)- 
split interstitial, rather than the (llO)-split configuration 
found in Fe Interestingly, the first principles calcu- 
lations also revealed that the (lll)-oriented SIA migra- 
tion energy is extraordinarily small (< O.OleV), which 
explains the experimental observation of diffusion down 
to 4 K SIA migration in Fe and V must therefore 

differ in their microscopic mechanisms. 

We perform a series of molecular dynamics (MD) sim- 
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ulation of SIA migration in V using an improved inter- 
atomic potential for V [l^ (refit to experimental and first 
principles data [l3| to reproduce the stable interstitial 
configuration) to address this discrepancy. In particular, 
we examine SIA diffusion as a function of temperature 
to determine the SIA migration mechanisms. We find 
that while SIA migration in V is similar to that in bcc Fe 
in many respects, its temperature dependence is highly 
unusual, exhibiting strongly non-Arrhenius behavior and 
correlation effects. It is this non-Arrhenius behavior that 
is the focus of the present Letter. 

The SIA is introduced in the form of a stable (111)- 
split interstitial and equilibrated for 10 ps at fixed tem- 
perature using a Langevin thermostat. The simulation 
was then switched to a microcanonical (NVE) ensemble 
in order to study SIA migration dynamics. Simulations 
were run at temperatures between 100 and 2000 K in a 
cubic simulation cell of edge length 10 ao, where ao is 
the temperature dependent-lattice parameter (the linear 
thermal expansion coefficient was 8.4 x 10~^ K~^). The 
diffusivity was measured by averaging over several 1 ns 
simulations. The interstitial position was identified by di- 
viding the space into Wigner-Seitz (W-S) cells centered 
around each perfect crystal lattice site. Interstitials are 
located in W-S cells containing more than one atom. 

Representative trajectories of the (lll)-split intersti- 
tial center of mass, collected over the whole 1 ns simula- 
tions, are shown in Fig. ^ For each temperature, more 
than 1000 jumps were observed, where an SIA jump is 
the exchange of an atom between neighboring W-S cells. 
As seen in Fig.^ the interstitial migration mechanism is 
strongly temperature dependent. For low and intermedi- 
ate temperatures (100-600 K) the (lll)-split interstitial 
executes a fully one-dimensional (ID) random walk along 
a (lll)-direction during the 1 ns simulation, as shown in 
Fig. da). At T - 700 K, the (lll)-split interstitial be- 
gins to make infrequent rotations from one (111)- to an- 
other (lll)-direction. This results in a three-dimensional 
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FIG. 1: Typical trajectories of migrating SIAs for tempera- 
tures of (a) 300 K, (b) 700 K, (c) 900 K and (d) 1400 K. 
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FIG. 2: Frequency of rotations as a function of temperature. 
The solid line is an Arrhenius fit to the data with slope AEr = 
0.44 eV, preexponential factor uq — 1.3 x lO^^s"^. 



(3D) trajectory that consists of long ID random walk seg- 
ments with abrupt reorientations, as seen in Fig. H^b). 
As the temperature increases, the frequency of the rota- 
tion events increases and the lengths of the ID trajectory 
segments decrease. At high temperatures, the rotation 
events become very frequent, leading to nearly isotropic 
diffusion (Fig. Hid)). 

Although these trajectories appear to be qualitatitively 
similar to those re por ted for other bcc metals (i.e. Fe 
and Mo) millillTI (cf- Fig.U^d) and Fig. 5 in [H), 
they differ in the elementary migration mechanism. The 
stable interstitial is the (lll)-split configuration in V, 
but the (llO)-split configuration in Fe and Mo. In the 
Fe and Mo cases, the split interstitial sits in the (110)- 
orientation until it is thermally activated into one of the 
(lll)-directions where it can mig rate easily before re- 
turning to a (llO)-orientation [l^. There are no relax- 
ation events of this type in interstitial migration in V. 
Here, the stable (lll)-split interstitial migrates long dis- 
tances and only requires significant thermal activation to 
reorient or rotate. 

The temperature dependence of the rate of rotation 
of the split interstitials from one (111) direction to an- 
other, ujr, in V is shown in Fig. [3 The data is well 
described by a conventional Arrhenius fit of the form 
LUr = j/Q exp [— A£'r/A:BT], suggesting that rotation is a 
simple thermally activated process. The activation en- 
ergy, AEr = 0.44 eV, is consistent with the energy dif- 
ference between the (111) and (110) configurations com- 
puted in first-principles (0.35 cV) and static calculations 
using the new interatomic potential (0.4 eV) [l^ . 

The diffusivity D (solid symbols) of the (lll)-split 
interstitial is shown in Fig. |21 for a temperature range 
between 100 K and 2000 K. D was determined from 
D = {R^{t)) /2t for ID diffusion, where the mean squared 
displacement {R^{t)) was calculated from the trajectory 



using standard averaging procedures [13|. Note that 
although the rotations at higher T change the dimen- 
sionality of the diffusion path, they do not contribute 
to transport and hence the diffusion mechanism is al- 
ways one-dimensional. If the diffusivity were Arrhenius 
[D/al = uoe-x:p[-AEd/kBT]), the data in Fig. El would 
lie along a straight line. This is clearly not the case; Fig.El 
shows pronounced curvature - especially at high temper- 
ature. Although Arrhenius behavior is widely expected 
for diffusion in the solid state, it is clearly inapplicable 
here. 

Non- Arrhenius behavior can have several different ori- 
gins. The energy barrier could simply change with tem- 
perature as a result of thermal expansion, as has been 
argued for the self-diffusion in bcc metals via a vacancy 
mechanism 0, 0| . The magnitude of the observed de- 
viations from Arrhenius behavior is too large to attribute 
to such a mechanism. The existence of multiple reaction 
pathways with different energy barriers can also lead to 
curvature in Fig.|3 However, detailed examination of the 
atomic configuration during diffusion shows that there is 
no change in diffusion mechanism over the entire tem- 
perature range. Although rotations are first observed at 

700 K within the 1 ns duration of the simulations, 
strong deviations from Arrhenius behavior are observed 
already at lower temperatures. A third possibility is that 
the degree of correlations in the diffusion process (i.e., 
particle jumps retain memory and the random walk is 
non-Markovian) is temperature dependent. Indeed, ex- 
amination of the SIA trajectories show that the (111)- 
interstitial has a higher probability of jumping back in 
the direction from whence it came, rather than forward 
along the same direction. At high temperatures, by con- 
trast, this effect appears to be reversed. 

We quantified this observation by measuring a corre- 
lation factor for split interstitial diffusion /, defined as 



3 



T[K] 

2000 500 200 100 




20 40 60 80 100 120 

l/kBT[eV-i] 

FIG. 3: A plot of the diffusivity of the {lll)-split interstitial 
in the form suggested by the Arrhenius relation. The filled 
symbols (■) correspond to the measured diffusivity D and 
open symbols (o) to D normalized by the correlation factor 
/ (see text). The straight line is a low temperature fit to 
the D/f data, corresponding to AEd ~ 0.018 eV and a pre- 
exponential factor of //o = 1.5 x lO^^s^^ in the Arrhenius form. 
Statistical errors are of order the symbol size. The inset to 
the figure shows the variation of / as a function of l/ksT. 



f — 2D/Db, where Di, is the "bare" diffusion constant 
defined as Db = if^n. Here, n is the mean number of 
jumps/second and Iq = \/3ao/2 is the jump length (near- 
est neighbor distance in the bcc lattice). If the intersti- 
tial trajectory is described by a sequence of jump vec- 
tors li, the correlation factor / is alternatively given by 
/ = 1 -I- 2^"~^(/i • Iq)/Iq, i.e. / = 1 for an uncorrelated 
random walk. The inset to Fig. O shows the variation 
of / with l/keT. At low T (T < 600 K), SIA motion 
is anticorrelated (/ < 1). The effect of the correlations 
on the diffusivity can be isolated by plotting D/f rather 
than D in Fig. O (open symbols). The low temperature 
D/f data (T < 500-?^) lie along a nearly straight line 
with slope AEd = 0.018 eV. Hence, the temperature de- 
pendent correlation factor explains the relatively weak 
deviations from Arrhenius behavior at low T. The un- 
usually small migration energy is consistent with experi- 
mental observations and the first principle estimate 
[ll|. One possible origin of the anticorrelations is that 
the finite relaxation time of the local environment around 
the SIA becomes less important as the thermal energy of 
the system increases. 

As the temperature is increased beyond 300 K, the 
correlation factor rises quickly to a value greater than 
unity. / > 1 is very unusual and may be thought of as 
correlated interstitial hopping over several barriers with- 
out completely thermalizing in between. This interpreta- 
tion is consistent with the fact that the correlation cor- 
rected diffusivity, D/f, only yields a straight line at low 
T but not at high T. We note that where D/f is ris- 
ing quickly, the energy barrier AEd obtained in the low 
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FIG. 4: The diflFusivity D (■) from Fig. El as a linear func- 
tion of r, normalized by AEd- The thick solid line represents 
the diffusivity of a particle in a sinusoidal potential at finite 
temperature as described by Eq. 10 for a value of 7 = 0.1r~^ 
where r = ymo^/e. The dashed line is the free particle limit 
of the same model, D = kBT/emy. The inset shows an Ar- 
rhenius plot of the diffusivity of the particle in the sinusoidal 
potential as a function of e/ksT as a solid line and the straight 
dashed line has slope one. 

temperature regime (where the behavior is Arrhenius) is 
smaller than the thermal energy. This is the source of 
the multiple interstitial hops at high T. Conventional 
derivations of activated escape over barriers j^^] usually 
assume AEd ~> ksT . 

Rather than trying to apply the Arrhenius description 
to interstitial self-diffusion in this system, a more general 
model is the motion of a particle in a periodic potential at 
all temperatures (from ksT < AEd to fc^T > AEd). 
In the limit that the barrier height is completely negligi- 
ble relative to the energy of the heat bath (i.e., a free par- 
ticle), standard arguments predict D — kBT/rwy, where 
m is the mass of the particle and 7 a relaxation time scale 
associated with a velocity-dependent friction force. Since 
this free particle diffusion coefficient is expected to be a 
linear function of T, we replot the data from Fig. |31on a 
linear temperature scale (Fig.^J. In this representation, 
the data is nearly linear, albeit with weak curvature at 
low temperature. This suggests that SIA diffusion in V 
is free particle like at high temperature (D ~ T), but 
follows the normal Arrhenius, hopping dynamics at low 
T {D ^ e-^^o/ksTy rp^g deviation from linearity at 
low temperature and the deviation from Arrhenius be- 
havior at high temperature suggests that a cross-over is 
occuring between the free and hopping particle limits. 

In order to better understand this transition, we explic- 
itly consider a simple one-dimensional particle of mass m 
diffusing in a sinusoidal potential by numerically solving 
the Langevin equation 

mx — ~ e/2cos[x/(j] + r], (1) 

where 77 is a Brownian white noise. Inserting values of 
m and a for V and e = AEd, this model yields nearly 
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linear difFusivity for 1 < ksT/e < 10 for reasonable val- 
ues of 7, followed by a crossover into Arrhenius behav- 
ior (see inset of Fig. 0)) for ksT/e < 1. Changing the 
value of 7 shifts the temperature at which the transi- 
tion from the particle hopping to free particle behavior 
is observed. This excellent agreement between the MD 
results and model predictions demonstrates that the ob- 
served strongly non- Arrhenius interstitial diffusion in V 
is a direct consequence of the relative magnitudes of the 
activation energy and the thermal energy. 

It is interesting to compare the situation described here 
for V to that for bcc Fe. The crowdion mechanism en- 
abling the easy interstitial migration along (111) direc- 
tions in V is also available in Fe, and estimates for the 
migration barrier AEd = 0.04 eV are similar to the V 
case. However, measurements of the apparent activation 
energy for diffusivity, analogous to the one presented here 
(albeit over a smaller temperature range between 700 K 
and 1200 K) yield much lar ger values of AEd = 0.12 
eV [iil or AEd = 0.17 eV 0|. This larger effective 
barrier is due to the fact that the (lll)-split interstitial 
must be thermally excited from the (llO)-state, i.e. the 
easy-diffusion configuration is not populated at all times 
as in V. The fraction of time during which the (111)- 
split interstitial is available for transport is given by = 
P{iii)/{P{iii) +P{iw)), where P^m) ~ exp [-AEf/keT] 
and P{iio) ^ cxp [— Ai^b/fesT] are probabilities for the 
interstitial to transform from the (110) state to the (111) 
state and back, respectively. Therefore, the interstitial 
diffusivity in Fe is the product of the diffusivity that the 
interstitial would have if it was always in the (111) ori- 
entation (hke in V) and F. For AEf > AE^, the effec- 
tive activation energy for interstitial diffusivity in Fe is 
~ [AEd + AEf), but for AEb > AEf, it is - AEd- 
Clearly, this imphes for Fe that AEf '3> AEf,. 

Molecular dynamics simulations of self-interstitial dif- 
fusion in bcc V were performed over an unusually 
wide temperature range (100-2000 K). Interstitial atoms 
in the (lll)-split configuration migrate very fast one- 
dimcnsionally along (111) directions during the 1 ns sim- 
ulations. As T is increased above 600 K, rotations of the 
split-interstitial from one (111) orientation to another 
occur with increasing regularity. The rotations can be 
described by Arrhenius kinetics with activation energy 
AEr = O.MeV. At temperatures T < 600 K, the diffu- 
sion exhibits significant anticorrelations. An Arrhenius 
analysis of the data (corrected for these anticorrelations) 
yields a very small migration energy barrier AEd = 0.018 
eV. For T > 600K, AEd is much smaller than the ther- 
mal energy and the Arrhenius expression is no longer 
applicable. The diffusivity then crosses over from Ar- 
rhenius to free particle type diffusion with increasing T. 
The fact that this type of diffusion is a linear function 



of T rather than Arrhenius, as usually assumed, could 
have important implications for predicting the lifetimes 
of reactor components in vanadium and other bcc metals. 
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